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A search for mixing in the neutral D meson system has been performed using semileptonic D° — > 
e + u decays. Neutral D mesons from D* + — > D°n + decays are used; the flavor at production is 
tagged by the charge of the slow pion. The measurement is performed using 253 fb~ 1 of data recorded 
by the Belle detector. From the yield of right-sign and wrong-sign decays arising from non-mixed 
and mixed events, respectively, we estimate the upper limit of the time-integrated mixing rate to 
be r D < 1.0 x 10~ 3 at 90% C.L. 

PACS numbers: 14.40.Lb,13.20.Fc,12.15.Ff,11.30.Er 



I. INTRODUCTION 

The phenomenon of mixing has been observed in the 
K" -K° and B° -B° systems, but not yet in the D° -~D° 
system. The parameters used to characterize D° — D° 
mixing are x — Am/T and y = Ar/2r, where Am and 
Ar are the differences in mass and decay width between 
the two neutral charmed meson mass eigenstates, and T 
is the mean decay width. The mixing rate within the 
Standard Model is expected to be small the largest 
predicted values, including the impact of long distance 
dynamics, are of order x ^ y ~ 10~ 3 — 10~ 2 . Observation 
of a mixing rate significantly larger than predicted would 
indicate either new physics (enhanced x) or insufficient 
understanding of long distance effects (larger y). 

For i,^ < 1 and negligible CP violation, the time- 
dependent mixing probability for semileptonic D decays 
isl 

V{D° -^D°^ X+tTVi) oc r D t 2 er Tt , (1) 

where rrj is the ratio of the time- integrated mixing to the 
time-integrated non-mixing probability: 

J~dtV(D Q ^D°^X+£-V t ) x 2 +z/ 
rD J™dtV(D°-,X-£+v e ) ~ 2 ■ 1 ' 

In this paper we present a search for D° — D° mix- 
ing using semileptonic decays of charmed mesons. The 
flavor of a neutral D meson at production is determined 
from the charge of the accompanying slow pion (jr s ) in 
the D* + — > -D°7r+ decay Q- We reconstruct the D° as 
D° — > K^~e + v e and identify its flavor at the time of de- 
cay from the charge correlation of the kaon and the elec- 
tron. We search for mixing by reconstructing the "wrong 



sign" (WS) decay chain, D*+ -> D°tt+, D° -> D°, 
D° — > K^ + e^V e , which results in a WS charge com- 
bination of the three particles used to reconstruct the 
candidate. The non-mixed process results in a "right 
sign" (RS) charge combination, ir+K~e + . In contrast to 
hadronic decays, the WS charge combinations can occur 
only through mixing, and t\d can be obtained directly as 
the ratio of WS to RS signal events. 

We make no attempt to reconstruct K*~ mesons. We 
treat charged daughter particles (K~ or 7r~) as if they 
were the direct daughters of the D°, omitting the ac- 
companying neutral particle. A small contribution from 
D° — » p~e + v e and D° — * ir~e + v e decays, due to charged 
pions surviving the kaon selection, is included in the anal- 
ysis. 

We measure rrj in a 253 fb -1 data sample recorded 
by the Belle detector at the KEKB asymmetric-energy 
e + e~ collider pj, at a center-of-mass (cms) energy of 
about 10.6 GeV. The Belle detector |5j] is a large-solid- 
angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov counters 
(ACC), a barrel- like arrangement of time-of- flight scintil- 
lation counters (TOF), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect K*l mesons and to identify 
muons (KLM). Two different inner detector configura- 
tions were used. The first 140 fb -1 of data were taken 
using a 2.0 cm radius beampipe and a 3-layer silicon ver- 
tex detector, and the subsequent 113 fb -1 were taken us- 
ing a 1.5 cm radius beampipe, a 4-layer silicon detector 
and a small-cell inner drift chamber 61. 



3 



Simulated events are generated by the QQ generator 
and processed with a full simulation of the Belle detector, 
using the GEANT package 0. To simulate mixed D 
meson decays we use generic (non-mixed) Monte Carlo 
(MC) events and appropriately reweight the proper decay 
time distribution. 



II. SELECTION AND RECONSTRUCTION 

For 7r+ candidates we consider tracks with: momen- 
tum p < 600MeV/c; projections of the impact parame- 
ter with respect to the interaction point in the radial and 
beam directions dr < 1cm and \dz\ < 2 cm, respectively; 
and electron identification likelihood ratio, based on the 
information from the CDC, ACC and ECL |j, C e < 0.1. 
Electron candidates are required to have p > 600MeV/c 
and C e > 0.95. Kaon candidates are chosen from the 
remaining tracks in the event with p >800MeV/c and 
C{K±) / (C^) + C^)) > 0.5, where £(£±,7^) is the 
likelihood that a given track is a or a 7r based on 
information from the TOF, CDC and ACC. Charged pi- 
ons originating from semileptonic D a — > ir^e + h' e decays 
that pass the kaon selection criteria are treated as kaons 
and assigned the kaon mass. According to MC simula- 
tion, this selection retains about 58% (43%, 52%) of sig- 
nal slow pions (electrons, kaons), and at this stage about 
13% (5%, 44%) of the selected ir s (e, K) candidates are 
misidentified. 

The variable used to isolate the signal events is 
AM = M(jt s Kev e ) - M(Kev e ), where M{n s Kev e ) and 
M(Kev e ) are the invariant masses of the selected charged 
particles and reconstructed neutrino (see below), with 
and without the slow pion. 

We require the cms momentum of the kaon-electron 
system to be p C ms(Ke) > 2 GeV/c to reduce the combina- 
torial background and the background from BB events, 
and to improve the resolution in AM. To further sup- 
press the contribution of BB events, we require the ratio 
of the second to the zeroth Fox- Wolfram moments, known 
as i? 2 , to be greater than 0.2 0- 

Background from D° — » K~tt + decay is suppressed 
by requiring the invariant mass of the K-e system to 
satisfy M(Ke) < 1.82GeV/c 2 . Background from D° -> 
K~K + decay is effectively reduced by the requirement 
\M(KK)-m D o \ > 10MeV/c 2 , where M(KK) is the K- 
e invariant mass, calculated using the kaon mass for both 
tracks. 

The background from photon conversions, where the 
e~ and e + are taken to be the electron and slow 
pion candidates, is suppressed by requiring M(e + e~) > 
150MeV/c 2 , where M(e + e~) is the invariant mass of 
the 7r+-e~ system with the electron mass assigned to 
the pion. A further search for an additional with 
charge opposite to that of the electron (slow pion) candi- 
date in the D° (D* + ) decay is performed in each event. 



Again, candidates with M(e + e~) below 150MeV/c 2 are 
rejected. 

The first approximation to the neutrino four-momen- 
tum is calculated from four-momentum conservation: 

Pu ^cms P-ir^Ke -Pest 7 (3) 

where P cms denotes the cms four-momentum of the e + e~ 
system, P^Ke the cms four-momentum of the ir s -K-e 
combination, and Pest the cms four-momenta of all re- 
maining charged tracks and photons in the event. The 
AM resolution is significantly improved by applying two 
kinematic constraints to correct P res t- First, events with 
-4GeV 2 /c 4 < M{ir s Kev) 2 < 36GeV 2 /c 4 are selected 
and the magnitude of P res t is rescaled by a factor x such 
that M{-n s Kev) 2 = (P cm s -jcPest) 2 = m 2 D , ± , with m D ,± 
fixed to its nominal value jlfj| . Next, only events with 
-2GeV 2 /c 4 < Ml < 0.5GeV 2 /c 4 are retained and the 
direction of the three-momentum p lcst is corrected in or- 
der to yield M 2 = 0. The correction rotates p rest in the 
plane determined by the vectors p Ics t and p^Ke- 

The mass difference AM is then calculated using this 
reconstructed neutrino four-momentum. The full width 
at half maximum (FWHM) of the AM peak in MC sim- 
ulated D° — > K~e + v events is found to be 7MeV/c 2 . 
Events with AM < 0.18GeV/c 2 are retained for further 
analysis. 

According to MC simulation, the D° — > K*~e + v e , 
p~ e + v ei and ir~ e + v e modes have AM distributions sim- 
ilar to that from D° — » K~e + v e decay, but with a larger 
FWHM. The charge correlation is preserved if the slow 
pion decays into a muon, and the muon is taken as the 
slow pion candidate. These decays will be referred to 
as "associated signal" decays. Their fraction in the re- 
constructed MC signal sample is 14.9 ± 0.1%, with the 
quoted error being due to MC statistics only. 

According to MC, the backgrounds for both RS and 
WS decays consist mainly of combinatorial background 
(98.4% and 98.0%, respectively). In order to minimize 
systematic uncertainties, we model the shape of this 
background using data, following a method validated 
with MC events. We use candidates with the unphys- 
ical sign combinations 7r+ K~e~ and n^K + e + to model 
the combinatorial background in the RS sample. In the 
WS sample this background is mainly due to random 
7r+ candidates combined with the selected kaon and elec- 
tron candidates, and is modelled by combining and 
(K~ + e + ) candidates from different events. The corre- 
lated background (1.6% in the RS, 2.0% in the WS sam- 
ple) arises from a true 7r^" from a D* + decay combined 
with true or misidentified K and e candidates, both (in 
the RS sample) or at least one of them (WS sample) 
from the corresponding D° decay. This background is 
described using MC. 

FigureQ]shows the AM distribution for RS events. We 
perform a binned maximum likelihood fit to the distri- 
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FIG. 1: AM distribution for RS events. The data (back- 
ground) is represented by points (hatched histogram), and 
the result of the binned maximum likelihood fit by the solid 
line. 



bution, maximizing 
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(4) 



where Ni is the number of entries in the i-th bin and 
/i(AMj) is the expected number of events in this bin, 
given by 

MAM,) = Af[f s P s (AMi) + (1 - / s )n(AM«)]. (5) 

P 8 (Pb) is the signal (background) AM distribution ob- 
tained from MC (as described above); f s and J\f are 
the signal fraction and the overall normalization, respec- 
tively, and are free parameters in the fit. The hatched 
histogram in Fig.Q] shows the fitted background contri- 
bution. The fitted signal fraction is f s = 73.4±0.1% and 
the number of RS signal events is = 229 452 ± 597. 

We increase the sensitivity to D° — D° mixing by ex- 
ploiting the measurement of the D° proper decay time. 
The decay time is evaluated using the measured momen- 
tum of the meson and the distance from the e + e~ inter- 
action point to the reconstructed D° decay K-e vertex. 
Due to the shape of the KEKB accelerator interaction 
region, which is narrowest in the vertical (y) direction, 
the dimensionless proper decay time t y is calculated as 



m D o j/vtx - Hip 



CT D 



Pv 



(6) 



where p y is the y component of the D° candidate's mo- 
mentum and y v tx and yip are the y coordinates of the 
reconstructed K-e vertex and of the interaction point, 
respectively. m D o and t^o are the nominal mass and 
lifetime of D° mesons [To| . 



Mixed and non-mixed processes have different proper 
decay time distributions: i 2 e~ t / TD ° and e - '/ TD ° , respec- 
tively. In order to calculate the mixing parameter 
after any selection based on D° proper decay time, the 
ratio of the two types of signal events, -/Vws/Wrs, must 
be multiplied by the ratio of their efficiencies, crs/ews- 
The efficiencies are obtained by integrating the convolu- 
tion of the above proper decay time distributions with 
the detector resolution function, over the selected t y in- 
terval. As the observed proper decay time distribution 
is a convolution of the signal and background probabil- 
ity density functions (p.d.f.) with the detector resolution 
function, the latter is found by performing a binned \ 2 
fit to the RS event t y distribution: J °° dt \ i f s e~ t / T D < > + 
(1 - / s )(/ e e-*/^ + (1 - f e )S(t))] x TZ(t/r D o - t y ). The 
signal fraction f s is obtained from the fit to the AM dis- 
tribution. The fraction f e of the non-prompt background 
component and its lifetime Tbkg are fixed to the values 
obtained by fitting the MC background t y distributions; 
S(t) describes the shape of the prompt background com- 
ponent. The resolution function lZ(t — t y ) is described 
phenomenologically by the sum of three Gaussians and 
an additional term for badly reconstructed tracks ( "out- 
liers"); we obtain the widths and coefficients from the 
fit. 

The ratios €rs / ews are given in Tabled The errors are 
obtained by varying each parameter in the proper decay 
time fit by ±lcr, repeating the fit and recalculating the 
ratios; the resulting changes are summed in quadrature. 
Using this method, the majority of systematic errors due 
to the imperfect description of the decay time distribu- 
tion cancel out. 

Comparison of the expected t y distribution for WS sig- 
nal and background events indicates that the figure-of- 

merit, defined as N^ s / J , is optimal for t y > 1.0. 
Events that satisfy this condition are retained for further 
analysis. 

III. FIT AND RESULTS 

We extract the RS and WS signal yields separately in 
six intervals of t y . The yields are obtained from binned 
maximum likelihood fits to AM, described by Eq. the 
WS signal distribution is taken to be the same as the RS. 
The results are shown in Table I and Fig. We obtain 
the time-integrated mixing probability ratio in the i-th t y 
interval, 



"d = AW^Vrs 



x e Rs/ e ws> b y multiplying the 
ratio of WS to RS signal events in each interval by the 
t y efficiency ratio. The mixing probability is compatible 
with zero in all proper decay time intervals (see Fig.[3J). 

The overall follows from a \ 2 fit with a constant to 
the measured r l D values, using a constant: 



r D = (0.20 ±0.47) x 10" 



(7) 



■5 




0.14 0.16 0.18 0.14 0.16 0.18 

AM [GeV/c 2 ] AM [GeV/c 2 ] 



TABLE I: The number of fitted signal events in the RS and 
WS samples, the efficiency ratio eRs/ e wSi an d the resulting 
r l D value for each proper decay time interval. 

t v ^Vrs ^ws 4sAws rh [1(T 3 ] 

1.00-1.33 18 742±166 -63.7±30.2 1.62±0.11 -5.49±2.63 

1.33-1.67 15 032±147 40.3±29.9 1.14±0.08 3.05±2.27 

1.67-2.17 16 430±155 -1.3±30.6 0.79±0.05 -0.06±1.48 

2.17-3.00 16 691±157 17.0±33.1 0.51±0.03 0.51±1.00 

3.00-5.00 15 443±160 14.7±37.4 0.30±0.01 0.28±0.72 

5.00-10.0 8 263±123 3.0±35.2 0.28±0.02 0.10±1.19 
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FIG. 3: Measured r l D values (points with error bars) in six 
different proper decay time intervals. The solid (dashed) line 
shows the null value (the fit to the six r l D using a constant). 
The dotted lines denote the statistical error of the fit. 



FIG. 2: WS AM distributions in six proper decay time in- 
tervals for data (points with error bars) and the results of the 
fit described in the text (solid line). At the bottom of each 
figure the fitted (expected for r % D = 90% C.L. upper limit) 
signal yield is plotted as a solid (dashed) line. 



The quoted error is statistical only. 



IV. SYSTEMATIC UNCERTAINTIES 

The experimental procedure was checked using a dedi- 
cated MC sample of mixed D° decays. These were added 
in different proportions to the generic MC, which includes 
both non-mixed D° decays and all known types of back- 
ground decays. The application of the same method used 
for the data on these samples verified that the recon- 
structed m value reproduces the input value without any 
significant bias. 

The main source of systematic error is the limited 
statistics of the fitting distributions, predominantly the 
background distribution in the WS sample. To estimate 



this uncertainty, we vary the contents of all bins of the 
RS and WS P S (AM) and P b (AM) distributions inde- 
pendently in accordance with each bin's statistical uncer- 
tainty, repeat the fit to the RS and WS data, calculate 
the corresponding r l D in each proper decay time interval, 
and obtain a new rjj value. Repeating the procedure, the 
obtained distribution of rjj values has a Gaussian shape 
with a width of 0.12 x 10~ 3 , which is taken as the sys- 
tematic error due to the limited statistics of the fitting 
distributions. 

To estimate the systematic uncertainty due to the AM 
binning and the shape of the WS background, the AM 
distributions in each proper decay time interval, origi- 
nally divided into 45 bins in the range 0.135 GcV/c 2 < 
AM < 0.180 GeV/c 2 , are re-binned into 12, 15, 20, 30 
and 60 bins in the same range. The fitting procedure is 
repeated for each set of bins, and the systematic error 
due to binning is taken as half the difference between the 
largest and the smallest m- ±0.07 x 10~ 3 . 

According to MC, the majority of the correlated WS 
background arises from D° — > K^tt + tt° decays. We 
repeat the WS fits, changing the amount of this back- 
ground in accordance with the relative error on B(D° — > 
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K~TT + ir°). The resulting variation in m is ±0.02 x 10~ 3 . 
The same method applied in the fit to the RS sample in- 
troduces a negligible variation of the rjj value (less than 
0.001 x 10~ 3 ). 

We compare the e^g values with the corresponding 
N^s/Njgg ratios. The relative difference between the two 
is conservatively assigned as the uncertainty on e^g/e^yg, 
although one expects at least part of the systematic un- 
certainty to cancel in the efficiency ratio. We reduce the 
six efficiency ratios simultaneously by this uncertainty 
and repeat the ro calculation; we then increase the ra- 
tios by this uncertainty and again recalculate rjj. We 
quote the difference between the resulting rjj values and 
the default fit (±0.02 x 10~ 3 ) as the systematic error 
due to an imperfect description of the proper decay time 
distributions. 

The systematic error due to the uncertainty in the as- 
sociated signal fraction is estimated by varying the frac- 
tion and repeating the fitting procedure. Using the errors 
on the measured branching fractions [HS of the associ- 
ated signal decay channels, we conservatively vary the 
amount of associated signal by ±40%. We recalculate 
rr> and compare it to the default ru value; we quote the 
difference, 0.002 x 10~ 3 , as the systematic error from this 
source. 

The quadratic sum of these individual contributions 
yields a total systematic error of ±0.14 x 10~ 3 . 

V. CONCLUSION 

In summary, we have searched for D° — D° mixing in 
semilcptonic D° decays. We independently measure the 
mixing ratios r l D in six D° proper decay time intervals. 
By fitting a constant to these six values we obtain the 
final result 

r D = (0.20 ±0.47 ±0.14) x 10~ 3 , (8) 

where the first error is the statistical and the second the 
systematic uncertainty. As we do not observe a signifi- 
cant number of WS D° decays, we obtain an upper limit 
for the time-integrated mixing rate. We use the Feldman- 
Cousins approach [llj to calculate the upper limit in the 
vicinity of the physics region boundary (td > 0). Using 
the result and the total error in Eq. [S] this yields 

r D < 1.0 x 10~ 3 at 90% C.L. (9) 



This result is the most stringent experimental limit on 
the time-integrated D° mixing rate obtained to date from 
semileptonic D° decays [To| . 
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